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Abstract 



£> ' We demonstrate that scattering of particles strongly interacting in three dimensions (3D) can be 

in ' 

suppressed at low energies in a quasi-one dimensional (ID) confinement. The underlying mechanism 
is the interference of the s- and p-wave scattering contributions with large s- and p-wave 3D 



scattering lengths being a necessary prerequisite. This low-dimensional quantum scattering effect 
might be useful in "interacting" quasi- ID ultracold atomic gases, guided atom interferometry and 
impurity scattering in strongly confined quantum wire-based electronic devices. 
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In a scattering process the boundary conditions play a key role and can substantially 
affect the scattering outcome. In case of asymptotically free motion and for spherically 
symmetric potentials, it is natural to analyze the scattering processes using partial waves 
belonging to certain angular momenta. Consider then a scattering event in 3D between two 
distinguishable particles with relative momentum k. Adding contributions of higher partial 
waves with increasing orbital angular momentum I and phase shifts Si can only increase the 
total cross-section, which in 3D is given by 
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(i) 



A pure s-wave (I = 0) approximation may not be accurate even at low energies for potentials 
that possess a significant spectral structure, e.g., impurity or dopant scattering in some 
bulk solid state systems or p-wave (1 = 1) interactions in ultracold atomic physics 

3, H Q , E , y, [3| ■ Such additional scattering terms for different partial waves could lead 
to increased heat dissipation along bulk conduction channels or to higher phase degradation 



HQ 



in atom interferometry I 111 . Il2t Il3l due to atomic collisions. 

The scattering process can change dramatically if it takes place in a partially confined 
geometry . su ch as quasi- ID systems. The latter situation is encountered in certain quantum 



wires 



u 



Of particular interest is the single mode regime, where only the ground state 
(hereafter indicated by the index 0) with respect to the confining directions can be actually 
populated. Then, instead of a, the quasi-lD total cross-section is the reflection coefficient 
R, which can be written in terms of the effective quasi-lD forward (z — > +oo) scattering 
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for electronic transport or in wave guides for ultracold atoms 



amplitude f$ 
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R = l-\l + f+ 



fo — fog + /( 
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(2) 



where fo g and fo u are even and odd quasi-lD amplitudes containing even and odd partial l- 
waves, respectively. In a seminal work |l8j], in the s-wave approximation for which f 0u « 0, a 
resonance R ~ 1 has been predicted with fo g w — 1 . The corresponding low energy transport 
is thus almost blocked and the longitudinal motion exhibits a near infinite effective quasi- 
lD scattering: the original confinement induced resonance (CIR) Q]. An account of the 
center of mass motion of colliding particles in the pure s-wave zero-range approximation is 
provided in 23]. On the other hand, higher partial wave contributions like the next p-wave 



2(j l22l ] in f 0u can also exhibit a resonance |20| with f 0u fa — 1 . This is another 
diverging effective interaction, but in the odd sector [soj ] . 

The present work addresses the novel situation where a simultaneous resonance in the s- 
and p-wave scattering occurs such that R fa 0. We show that the even and odd contributions 
to the scattering amplitude interfere and may provide effectively a quasi-lD gas of non- 
interacting atoms or a free flow of carriers in a quantum wire. Note, that this free flow 
occurs in spite of a strong interaction in 3D. We remark that for cold atoms dc electric field 
can supply us with hybridized Feshbach resonances providing a situation where both s- and p- 
wave scattering becomes simultaneously resonant [2jJ . The above-mentioned complementary 
or "dual" CIR due to both large s- and p-wave scattering is established by extracting the 
transmission T = 1 — R = \1 + f^\ 2 from the wave packet dynamical solution of the time- 
dependent Schrodinger equation for two particles with masses mi ^ m 2 and by highlighting 
the key concepts at play. The external transverse confinement is assumed to be parabolic 
Ui(pi) = \miOjfpf where 7*j = (pi,Zi), i = 1,2, are the coordinates of the particles. Their 
interaction is the screened Coulomb potential (Vq < 0) 

V(r) = V je- r ^, (3) 

where r = |t*i — r%\ and tq is the screening length. Note that only for uj\ = u 2 the center of 
mass decouples and e.g. an impurity-like scattering can be described solely by the relative 
coordinates. If lj 1 ^ lu 2 , a two-species gas of ultracold trapped atoms, e.g., 40 K, 85 Rb, 
87 Rb and combinations thereof, can be described qualitatively, since only the low-energy 
asymptotic behaviour of the scattering process is required, which should depend mostly on 
properties such as the scattering lengths and the range of the potential V. Independent 

n n 

simulations with other potentials, e.g., Cq-C\ 2 PJ| or hyperbolic-cosine |2£J] are expected to 
provide the same behaviour and properties. Additionally Eq.fjHl) allows for a good control 
of the low energy scattering lengths by tuning its spectrum via Vq. The computation of T 
is first performed for U\ = u 2 and second for ui ^ u 2 , including the coupling to the center 
of mass. 

In the Schrodinger equation, the center of mass coordinates R = + ^-r 2 , with 

M = mi + m 2 , are employed in the cylindrical representation (pa, (pa, Z) while for the 
relative coordinates r = T\ — r 2 , depending on the representation both cylindrical (p,<f),z) 
and spherical coordinates (r, 6, <fi) are needed. By symmetry, the motion with respect to the 



degree of freedom Z can be decoupled from the problem and will be omitted. In addition, due 
to the conservation of the ^-component of the total angular momentum, namely, L l2 + L 22 = 
—%{dfa + dfc) = —i(d^ + d<j, R ), (pi being the azimuthal angle of Tj (h = 1), the number of 



mass 



independent variables can be further reduced in a frame co-rotating with the center o 
around the symmetry z-axis. This is accomplished by the unitary transformation 

W = e^ iz , (4) 

where L z = —id^. Expressing the total Hamiltonian in terms of pR,<j)R,r,p,9 and <p, the 
net effect of this transformation is to shift d ( f )R to d$ R — in the kinetic operator for the 
center of mass and (p to (p + (pR in the potential operator that couples p to pr. The rotated 
Hamiltonian without V becomes then 

H = H M + H fl + W. (5) 

Here H M = -(d 2 R + l/Ap 2 R )/2M - (d^ - d^/2Mp R , W = \Mu\ lP \ + \p^p 2 + p(u 2 - 
uj^ppRCOscp and = —d 2 /2p + L 2 /2pr 2 , where L 2 is the square of the orbital angular 
momentum of the relative coordinates with eigenvalues 1(1 + 1), p = mim^/M is the reduced 
mass, u 2 = (m2/M)u}f + (mi/M)uj% and u 2 M = (m 1 /M)(J 2 + (m 2 /M)uj% characterize the con- 
finement oscillator frequencies for the relative and center of mass coordinates, respectively, 
whereas V is invariant with respect to the transformation U. The derivative in the sec- 
ond term of H M is a Coriolis coupling term, which is characteristic for the rotating frame of 
reference and appears as an additional kinetic energy operator. 

The initial wave packet at t — possesses a Gaussian shape of width a z and momentum 
ko = \J2pe > for the unconfined ^-motion 

(0) = Nr^ e - {z - zo)2/2a * e izk °, (6) 



N is the overall normalization constant defined as (\E'(0)| v I / (0)) = 
J °° dr J °° dpR J* sinOdO dcp^* (0)\I/(0) = 1 and the particles are separated by z < 
while their confined motions pi are in the respective harmonic oscillator ground state 
$ = $(pi,p2) = e~( p i/ a i +p 2/ a 2)/ 2 with di = (1 1 miUJi) 1 ! 2 . By expressing pi in terms of the 
variables pr, 4>r, p and 0, the action of U is to simply shift <p to (p + (pR in the exponent 
of $, as can be seen by using the basis of eigenstates of L z . As a result, the cpR variable 
drops from the initial state \]/ (0) and the Schrodinger equation i^>(t) = [H + V] ^(t) 



4 



reduces to the subspace of the four variables (pu, r). The z— component of the total angular 
momentum is in the co-rotating frame given by Lz = —id^ R . According to our initial 
wave packet in Eq.© we assume now a vanishing total angular momentum Lz^ = and 
therefore omit the corresponding derivative terms in Hm- 

The extraction of T from follows from the asymptotics at large times t, when ^(O) 
has been split by V into a backward and a forward scattered part outside the range of V 
ideally for z — ► =Foo, respectively. It can be obtained from the forward scattered part by 
integrating the asymptotic density |^/(t)| 2 over the half-space z > r or from the projection 
|(\I/o(^)| x I / (^))| 2 — ► |1 + fo\ 2 onto the unscattered large t evolution of ^(0) under H alone. 

t— >oo 

In the limit a z — » oo our wave packet scattering corresponds asymptotically (t — > oo) to the 
monoenergetic stationary scattering situation thereby providing the energy or momentum 
resolved scattering parameter T = T(k). 

The present computational scheme to solve for is a four dimensional extension of 
a three dimensional method originally developed to treat, amon g ot hers, bound-bound and 
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28| . In this scheme, 



bound-continuum transitions for atomic systems in external fields 
an angular basis fj(9,<p) is constructed on the grid (9j,4>j) using the exponentials e im ^ and 
the Legendre function P/™(#), such that the only non-diagonal terms are the angular parts 
of the kinetic energy operators, namely, L 2 /2/zr 2 and — <9 2 /2Mp|,. These, in turn, can be 
diagonalized efficiently if one employs additionally a simple unitary transformation between 
the two basis fj and e im<?i , P^iO) I^. Our approach is reminiscent of a two dimensional 
discrete variable representation (DVR) [29] and the values of the wave function ^(t) are given 
on an angular grid (9j,<j)j) whose total number of grid points equals the number of basis 
functions fj. The solution is then propagated in time using a component-by-component split- 



operator method 
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28j according to (t + At) w e -iwAt e -i(H^+v)At e ~iH M At e -^wAt^,^_ 



To be specific the parameters are chosen as follows. r = 1 defines the length scale, 
mi/ m2 = 40/87 is the mass ratio (e.g., 40 K and 87 Rb) and the reduced mass is taken as the 
unit mass /J, — 1. We define u = (u>i +002)/^ and the confinement length a± = sJYfyio. 
For the case of a decoupled center of mass motion, uji = uj 2 = uj, whereas for the coupled 
case, LUi = 1.35u;2. The low longitudinal energy e = 0.002 guarantees the single mode regime 
u<E = uj + e<3uj for the region of 0.002 < uo < 0.02 we consider. 

For — Vq < 9 (in units of h 2 / T^o), there exists no bound-state of V(r) + 1(1 + l)/2/ir 2 for 
I > 2. Thus only the s- and p-wave 3D scattering lengths a s and a p can be large. These 
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are computed in the limit k — > by kcot5 (k) = —l/a s and k 3 cot5i(k) = —l/V p , with 

V P = 4. 

First let us consider the confined scattering process for the case of a decoupled center of 
mass motion uo\ = uo 2 = u. The results for T are shown in Fig.^together with the behaviour 
of the 3D scattering lengths a s and a p . These are tuned by changing Vo in two regions where 
correspondingly only a s (—V ~ 1) or both a s and a p (—V ~ 7.5 — 9.5) are large , i.e., 
on the order of a±. In the first region we see the well-known s-wave CIR which leads to a 
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3- 



minimum T pa when the ratio a±/a s approaches a±/a s pa 1.46 in agreement with 
In the second region, however, remarkable peaks T pa 1 in the transmission are observed 
being most pronounced for small momenta ko (small a±). Considering in this case pure 
s-wave scattering would yield T pa in contrast to the behaviour T pa 1 shown in Fig. d 
Clearly this is due to the fact that both a s and a p become large and contribute equally to 
the scattering process. We therefore encounter the peculiar situation of an almost complete 
transmission, i.e. free flow, in spite of the strong interaction between the scattering partners 
in free space. 

This effect of the suppression of the quantum scattering (T pa 1) for —Vo ~ 8 — 9 is 
encountered equally in the case of a coupling with the center of mass U\ 7^ uj 2 (see Fig. |2J). 
The position a±/a s pa 1.45 of the total transmission for u}\ = uj 2 = 0.02 is only slightly shifted 
to a±/a s pa 1.25 under the action of the coupling with the center of mass for oj\ = 1.35co>2 
and the maximum of T is slightly decreased. Note that in strong contrast to T, the 3D cross 
section a, e.g. for Vq = —8.45 and E ~ uo\ = uo 2 = 0.02, can increase almost four-fold if the 
p-wave contribution is added. 

The observed dramatic change of T due to the combined action of s— and p-wave scat- 
tering is qualitatively confirmed in simpler though solvable models, at least in the absence 



of the coupling to the center of mass. Indeed, for u\ = oo 2 and kr = y/2jjEr 1 (i.e., 
low total energy), the condition equivalent to fo u pa —1 can be obtained separately from the 
antisymmetric part of the wave function as a divergence of a quasi-lD effective interaction 
strength for a±/a p = —1.36 J^l^. Together with a± = 1.46a s [18], one finds again the re- 
quirement of both large a, ^ Another more direet determination of both /„„ nsing the 
approximation developed in |22j predicts a± = 2a s = — 2a p , namely, f 0g = —(1 + icot5 g ) 1 
and fou = —(1 + icotc^)" 1 , where cot S g = — [a±/a s — (C 2 — a 2 L kQ) 1 ^ 2 ]a±k /2 while 
cotS u = —[a\/V p + (C 2 — a 2 ± kl) 3 ' 2 }/ (6a±k ) is a straightforward improvement to Eq.(30b) 
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FIG. 1: (color online) (a) The 3D scattering lengths a s (black) and a p (red) and (b) the quasi-lD 
transmission coefficients T for a few values of a± and fixed e = 0.002(&o = 0.0632). Discussion see 
text. 

of Ref. (22I and (7 = 2. Using then a s (V ) and a p (V ) of a square-well of depth V and radius 
r in the scattering amplitudes /o 9)U we obtain T from Eq.(J2J) which is shown in Fig. El The 
resulting behavior is qualitatively similar to Fig. ^ confirming the role of both large a s and 



"p 1 



on the order of a±. 



For an experiment, the main condition for T ~ 1 is /o s , u ~ — 1. If the partial waves 
I > 2 are negligible and fco r o kr <C 1, this results in the fact that a s and — a p are 



of the order of a_i_. More precise values of these ratios 



IS. 20] may depend on del ails 



of the confining potential, coupling to the CM, k , and should be computed for a given 
experimental setup. Current laser traps for e.g. 40 K can well reach the ranges a s ~ a± < 
60nm j^lj or — a p ~ a± < 60nm [ijj], usin g, h owever, separate 3D Feshbach resonances in 
the s- and p- waves p]. In a recent paper |24j a mechanism for the simultaneous creation 
of s- and p-wave Feshbach resonances is established: Applying laboratory dc-electric fields 
introduces the herefore necessary mixing of s- and p- waves. Thus, a reduction of quasi- 
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FIG. 2: The influence of the center of mass coupling on the transmission T for u>2 = 0.02 and 
e = 0.002. Here uj\ = 1.35(^2, a±Ao = 6.52 (solid curve) and oj\ = u>2, a±/ro = 7.07 (dotted 
curve) . 

ID scattering relative to 3D, as revealed e.g. in trap losses, should be observable with 
tighter traps by further exploring Feshbach resonances. Another example is the potential 
V = —(ZiC 2 1 ' nr)e~ r / r ° of an ionized impurity of valence Zi in a semiconductor quantum wire 
with dielectric constant k. Restoring the length unit 7"o, the relation to Vq (in units of h 2 / /ir^) 
is thus r = —Kfr^Vo/Zi/ie 2 . For Vq = —8.45, the required screening length r can range from 
5.4 nm (heavy holes in Si and Zi = 2) up to 86.8 nm (electrons in GaAs and Zj = 1). 
Geometrically, this range and a± 3> tq are well within current scaling technologies [12( . 

In conclusion, due to quantum interference of different partial wave amplitudes, the scat- 
tering in a quasi-lD geometry can be strongly suppressed, although the interaction or scat- 
tering in 3D is strong. 
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by P.S. and V.S.M. 
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